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The quality of seismic-to-well tie is commonly quantified using the classical Pearson's corre- 
lation coefficient. However the seismic wavelet is time- variant, well logging and upscaling is only 
approximate, and the correlation coefficient does not follow this nonlinear behavior. We introduce 
the Dynamic Time Warping (DTW) to automate the tying process, accounting for frequency and 
time variance. The Dynamic Time Warping method can follow the nonlinear behavior better than 
the commonly used correlation coefficient. Furthermore, the quality of the similarity value does not 
depend on the selected correlating window. We compare the developed method with the manual 
seismic-to-well tie in a benchmark case study. 
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Introduction 

Well logs are commonly used as ground truth to correlate the seismic signal with the earth's stratigraphy 
( White and Simirl 2003 ). In this process, a wavelet is first estimated from the seismic trace, and then it 
is convolved with the reflectivity calculated from the well logs (sonic log and bulk density log). Iterative 
techniques are used to estimate the wave let with correct phase and amplitude spectrum by matching the 
actual seismic trace at the well position ( Hampson-Russelll . 1999 ). Tying the seismic traces to the well 
logs aims to m inimize the differences in the way seismic data and well logs measure the same parameters 
( Burchl . 2002 ). but with different resolution. 



The quality of the tie between the synthetic and the seismic trace is based on the correlation coefficient, 
which is limited to linear features. The time- variant nature of the seismic wavelet adds nonlinearities to 
the trace which cannot be followed by a linear metric. Thus, wavelet phase mismatches frequently occur 
between the final processed seismic data and the synthetic seismograms created from well logs. Thi s 
fact leads to potential complications in stratigraphic and structural interpretation (I Van der Baanl . l2008r) . 
We propose a new method to match these time series accounting for frequency and time variance. We 
argue that Dynamic Time Warping (DTW) method can follow these changes and furthermore the quality 
of fit is not limited to the selected correlation window. 



Nonline arities in time series repr esenting physical processes are common in m any areas (speech pro- 
cessing (IRabiner and Juangl . 119931 ) , medicine, industry and fi nance dKeoghl.l2002h\ More recent l y, thes e 
concepts have been improved by iKeogh and Kasettyl (120031 ) and iKeogh and Ratanamahatana (12004). 
DTW is a robust tool to match time series even if they are out of phase or time shifted ( Keoghl . 2002b . 



Related wo rks using DTW in seismic applications co me with the attempts to automate well-to- well log 



correlation (ILineman et all 119871: IS teven et al., 2004). Well logs from different wells are correlated to 



infer common earth features. ISteven et alJ (120041 ) found that the cross-correlation was unable to follow 
local distortions such as st retching or shrinking of stra tigraphic intervals, typical of logs collected even 
from closely spaced wells. lAnderson and Gab 3 <[L983) seek for correspondence between features in the 



logs of various wells using dynamic programming tools. To our knowledge, there are no previous reports 
of dynamic programming applied to the seismic-to-well tie problem. 
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Method and Theory 



The c orrelation coefficient is commonly used to measure the quality of the seismic-to-well tie ( Hampson-Russelll . 



19991) . Comparing two (time-dependent) sequences S := [s\,S2, • • • ,s n ] and T := [t\,t2, . . . ,t n ], both of 



length n, will give a correlation coefficient at the time lag v. 
where ji x is the average of trace x. 

This measure works well if a constant time shift t characterizes both signal. B ut the majority of geophys 



ical applications have time alignment problems ( Anderson and Gaby . 1983b . When this time alignment 



is constant, the problem is reduced to the correction of the time lag by cross correlation. But this measure 
fails to find the best matching in nonlinear cases. 

An alt ernative to the cross correla tion is to find the Euclidean distance (L^— norm) between the two time 



series (IKeogh and KasettyU2003h : 



Deuciid(SJ) = 4/£(S(i)-r(i)) 2 - (2) 



i=\ 



where D euc nd(S, T) is the one-to-one distance between the synthetic S and the trace T. 

The Euclidean distance (L2— norm) is the most widely use d distance measure. It is tr i vial to implement 
but al so is very sensitive to small distortions in the time axis dKeogh and Kasetty . 2003 : Berndt and Clifford). 



994). Taking the advantages of the Euclidean distance and adapting it for nonlinear matching, Berndt and Clifford! 



(119941 ) proposed the Dynamic Time Warping technique as we know it today. 



DTW distance can accommodate stretching and squeezing in the time series by linear programming. It 
uses the Euclidean distance as the initial metric but allows for the one-to-many non-linear alignment. 
The warping distance is represented as the minimum path in a grid representation of both sequences. In 
Figure Q] the warping path, W = w\,W2, . . . ,w^, aligns the elements of S and T is such a way that the 
distance between them is minimized. 
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Figure 1 Warping path at the minimum distance of two sequences. The warped versions of this example 
are S warpe d = Sl,s2,s2,s3,s3,s4,s5,s6,s7,s7 and T warpec i = tl,t2,t3,t4,t5,t5,t5,t6,t7. 

In this matrix the square distance in the elements (i th J th ) is calculated by: 

8(s h t i ) = (s i -t i ) 2 (3) 
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To find the best alignment between these two sequences we have to retrie ve the path through the matrix 



that minimizes the total cumulative distance bet ween them (IKeoghl. 12002) as illustrated in Figure [TJ The 



optimal path minimizes the total warping cost (Be rndt and Clifford 1 19941) is: 



DTW(SJ) =min£S(4 (4) 



w 

w k=\ 



where each Wk corresponds to a point (ij)k- From Figured] we can extract the first samples of the new 
warped sequences as S warp ed = [^1^2,^2 ,£3 ^4] and T warpec i = [h^ ,^3,^3,^4]. In this way sequences are 
accelerated or decelerated along the time axis. From the linear programming point of view the problem 
is to find the minimum cost warping path, \i^ — h\ < Wfc The dynamic p rogramming approach uses the 



following recurrence to find the warping path dBerndt and Clifford! . 1 19941 ): 



y(i, j) = 8(si,tj) + min [y(i - 1 , j) , y(i - 1 , j - 1) , y(/, j - 1)] , (5) 

where 8{s^tj) is the distance defined in ©, and the cumulative distance y(/,j) is the sum distance 
between the current elements and the minimum cumulative distance of the three neighboring cells. 

We apply the DTW algorithm to obtain a first seismic-to-well tie alignment between observed seismic 
data and the synthetic trace created from the well logs. 

Results 



The dataset used in our experiments consist o f a 3D post- stack seismic profile, with 13 wells and their 
correspondent logs ( Hampson-Russelll . 1999b . We use well 08-08 to estimate the wavelet which is 



subsequently used in all ties. The original observed data and unstretched synthetics are subjected to the 
DTW approach. 

The estimated warping path for well 01-08 is shown in Figure [2] (left). Note the nonlinear relationship 
between both signals. 
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Figure 2 Left plot shows the warping path represented by the yellow nonlinear trace, the lateral and 
bottom signals are the Trace and the Synthetic respectively. The right upper plot shows the original 
signals and their corresponding matching points, note that nonlinear and one-to-many associations 
occur in both ways. The right bottom plot depicts the warped version of both signals. 

The effect of the DTW over the signals while finding the best match is shown in Figure [2] (right). 
Instead of having a one-to-one comparison between these two signals a nonlinear alignment between 
their matching points is found. By applying the time warping correction to the original signals we 
obtain their warped version as is shown in the lower right plot of Figure [2 They are both stretched to 
match their common features. 
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Figure 3 Case 1: Well 01-08. Upper plot shows the manual seismic-to-well tie, the correlation of this tie 
was 80 %. Bottom plot shows the automatic tie with a correlation of 83 %. A good agreement is found 
between both the manual and automatic well tie. 



The correlation coefficients are estimated over the entire length of the well log for the DTW. Results 
of this measure are similar to the ones obtained in the time window from 800 ms to 1100 ms for the 
manual tie. Figure [3] shows the results for the manual seismic-to-well tie for Case 1: well 01-08. The 
DTW method was able to match similar events along the full seismic trace and gives a similar correlation 
coefficient to the one obtained using the manual method. 

Figure |4| shows a second example where a better correlation is obtained in the automated method than 
the manual result. Note the correspondence between both traces in the bottom plot. 
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Figure 4 Case 2: Well 16-08. Upper plot shows the manual seismic-to-well, the correlation of this tie 
was 74 %. Bottom plot shows the automatic tie with a correlation of 80 %. 

Conclusions 

We have implemented dynamic time warping to automate the seismic-to-well tie procedure, as this 
approach provides an optimal solution for matching similar events. We strongly advocate however to 
beware fully automated and non-supervised applications of this method, as a visual quality control of 
the end result remains highly advisable. 

Future efforts will be oriented to analyze the effect of the estimated wavelet over the resultant warped 
sequences. Many other applications of DTW are envisionable for seismic data. These include log-to-log 
correlations, alignment of baseline and monitor surveys in 4D seismics, PP and PS wavefield registration 
for 3C data. 
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